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Braginsky, Gorodetsky, and Vyatchanin have shown that thermo-refractive fluctuations are an 
important source of noise in interferometric gravitational-wave detectors. In particular, the thermo- 
refractive noise in the GEO600 beamsplitter is expected to make a substantial contribution to the 
interferometer's total noise budget. Here we present a new computation of the GEO600 thermo- 
refractive noise which takes into account the beam's elliptical profile and, more importantly, the 
fact that the laser beam induces a standing electromagnetic wave in the beamsplitter. The use of 
updated parameters results in the overall reduction of the calculated noise amplitude by a factor 
of ~ 5 in the low-frequency part of the GEO600 band, compared to the previous estimates. We 
also find, by contrast with previous calculations, that thermo-refractive fluctuations result in white 
noise between 600 Hz and 39 MHz, at a level of 8.5 • 10 _24 Hz _1,/2 . Finally, we describe a new type 
of thermal noise, which we call the thermo-chemical noise. This is caused by a random motion 
of optically-active chemical impurities or structural defects in the direction along a steep intensity 
gradient of the standing wave. We discuss the potential relevance of the thermo-chemical noise for 
GEO600. 



INTRODUCTION AND MAIN RESULTS FOR 
THERMO-REFRACTIVE NOISE 



The optical layout of GEO600, the German-British 
gravitational-wave detector, differs in an essential way 
from that of LIGO and VIRGO, the two other work- 
ing gravitational-wave interferometers. In both LIGO 
and VIRGO each of the arms consists of a high-quality 
Fabry-Perot cavity, where most of the light power is con- 
centrated, while the beamsplitter is located outside of 
both cavities. By contrast, in GEO600 the light, after 
passing the beamsplitter, is processed only once through 
each of the doubly-folded arms pp. Therefore, GEO600 
is much more sensitive to the noise originating at the 
beamsplitter than its LIGO/ VIRGO counterparts. 

The beamsplitter in the GEO600 laser interferometer 
induces a difference in the optical path lengths between 
the interferometer's two arms since one of the beams is 
transmitted through the beamsplitter twice while is the 
other is being reflected from its surface. This difference 
in the optical path length can be modified by a change 
in the refractive index of the beamsplitter due to random 
fluctuations of temperature in the refracting material. As 
Braginsky, Gorodetsky, and Vyatchanin have shown [2J, 
the temperature fluctuations result from a random heat 
flow between different parts of the refractive medium. 
Since the interrogating beam is not homogeneous but 
has a Gaussian profile the inhomogeneous temperature 
fluctuations produce a net (fluctuating) change in the 
transmitted beam's phase which is seen as noise in the 
overall output of the instrument. For a circular Gaussian 
beam, Braginsky and Vyatchanin |3J have computed the 
following spectral density of this thermo-refractive noise 
expressed as an equivalent change in the end mirror's 



position 6z: 
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Here, ks is Boltzmann's constant, k the thermal con- 
ductivity, T the temperature, (3 = dn/dT where n is 
the refractive index, a is the beamsplitter's thickness, 
C the specific heat, p the density and r$ is the beam's 
radius defined in terms of the beam intensity, so that 
I = Iq cxp(— r 2 /r 2 ,). The expression in Eq. (JTJ) has 
been used to estimate the thermo-refractive noise in the 
GEO600 beamsplitter g]. 

However, inside the GEO600 beamsplitter, the beam is 
not circular: the 45° angle of incidence and the refractive 
index n = 1.45 of fused silica introduces ellipticity of the 
beam. The ratio of the major and minor axes of the 
beam's elliptical crosscction is given by T) = 1.23, while 
the minor axis equals the width of the original beam. 
More importantly, the beam is not homogeneous in the 
longitudinal direction since a standing wave is formed in 
the beamsplitter as a result of the beam's reflection from 
the end mirror. The changes in the material's refractive 
index [5] close to the nodes of the electric-field intensity 
have less effect on the wave's phase than changes in the 
antinodcs. Therefore, random thermal fluctuations on 
the scale of the light wavelength are important, making 
for a much higher level of noise at frequencies higher than 
~ 600 Hz than the BV result would suggest. 

In this paper we derive the expression for the thermo- 
refractive noise with these corrections taken into account. 
We obtain the following result: 
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Here a' = a/ cos(i) is the optical pathlength through the 
beamsplitter, where i — arcsin (1/V2n) ~ 29° is the angle 
between the beam axis and the normal to the beamsplit- 
ter, k = 27m/ A is the wavevector and Ith — \J k/(C pu) is 
the thermal diffusion length. The numerical values rele- 
vant for GEO600 are given in table [I] The first ter m on 



Symbol 


(GEO600-)Value 


k B 


1.38- 10" 23 J/K 


K 


1.38 W/mK 


T 


300 K 


P 


8.5- 10 -6 


a 


8.0 • 10~ 2 m 


V 


1.23 


k 


8.56 ■ 10 6 nf 1 


TO 


0.71 • 10~ 2 m 


c 


746 J/kgK 


p 


2200 kg/m 3 


L 


1200 m 



TABLE I: Parameter values taken from [B] 

the right-hand side of Eq. Q is the contribution due to 
the transverse elliptic profile of the beam, and it agrees 
with the result in [3] for r) — 1, while the second term is 
due to the presence of the standing wave. It is common to 
express the noise amplitude in terms of the dimensionless 
metric h for easy comparison to other sources. Figure [T] 
is a plot of VS^ — V S Sz /L where L is the interferome- 
ters unfolded armlength and figure [2] is a comparison to 
the measured noise spectrum and the previous estimate 
for the beamsplitter thermo-refractive noise [7J. 
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FIG. 1: loglogplot of the noise amplitude \/S h (f) in units of 
Hz -1,/2 as a function of / in units of Hz with the GEO-600 
parameters 

At a typical gravitational wave frequency of 100 Hz 
the noise amplitude is y/W* ~ 5.1 • 10~ 23 Hz -1 / 2 . For 
low frequencies the effect of the standing wave is negligi- 
ble and we have the 1/ / dependence, consistent with the 
previous calculations. However, for higher frequencies 



the standing-wave contribution takes over and the noise 
spectrum becomes white between 600Hz and 39Mhz with 
an amplitude VS^ ~ 8.5 ■ 10~ 24 Hz- 1/2 . At the latter fre- 
quency the thermal diffusion length becomes comparable 
to the wavelength of the beam light, and the 1// depen- 
dence is recovered but at much higher value than would 
be predicted by the BV formula. 

II. METHOD OF CALCULATION 

The foundation of this analysis is direct method [5] for 
calculating thermo-refractive noise. The relevant readout 
variable is the phase of the output beam, translated into 
a change in the end mirror's position Sz. It is a fairly 
straightforward calculation to show 
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where the integral is over the beamsplitter's volume V, 
5T(r) is the temperature fluctuation at r, and q{r) is the 
form-factor given, taking the elliptic Gaussian beam and 
standing wave into account, by 
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We now calculate the noise of this generalized vari- 
able via the hypothetical experiment proceeding in three 
steps: 

1. Periodically inject entropy into the medium with 
volume density 

SS(r, t) 
dV 



F cos(ojt)a'f3q(r) (5) 
as a result of 



Calculate the dissipated power Wd 
this entropy injection. 

To do the latter, we solve the heat equation 
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dt dt 
while keeping in mind that: 

• the solution is periodic with frequency uj 

• the wavelength is much smaller than the 
beamsplitter thickness so we can ignore 
boundary effects 

• the diffusion of the osscilating temperature in 
the transverse direction is negligible, giving 
V 2 = d 2 /dz 2 , since ro >> Uh for all frequen- 
cies of interest. 

The dissipated power is then computed by using 
the standard expression 
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taken from e.g. Landau and Lifshitz [5j, where 
denotes the time average over one period. 
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3. Calculate the spectral density of the noise in 5z 
using the fluctuation dissipation theorem [8], |10J : 
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These three steps lead directly to the main result in 
Eq. {2}. 



III. THERMO-CHEMICAL NOISE 

The small scale of the wavelength opens up the pos- 
sibility for a new type of thermal noise which has not 
yet been considered. The fused silica used for contem- 
porary beamsplitters contains minute quantities of con- 
taminants such as OH ions, CI ions and other defects 
which have an effect on the refractive index depending 
on their concentration. As these optically-active con- 
taminants diffuse up and down the steep gradient of the 
standing-wave electric-field intensity, they cause fluctua- 
tions of the overall beam's phaseshift. 

Lets compute this thermo-chemical noise. Let P(r) 
be the fluctuating volume concentration of the optically 
active impurities. The optical path change due to these 
impurities is given by 



Sz= / d 3 f- a' ■ a - 5P(f) -q(r), 
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where a = dn/dP and q(r) is given by Eq. [4] 

We can follow our earlier treatment of thermal noise 
[TT] and calculate the dissipated power in the system un- 
der the Hamiltonian 
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This can be easily done by recalling that the formal ex- 
pression for P(f) is 



P(f)=S i S('r-f i ), 



(11) 



where fi are the positions of individual optical impurities. 
Then, under the action of the above Hamiltonian, each 
impurity experiences a force 



t / ,\ / dq(fi) 

fi = F cos{ujt)a a — . 
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Under the action of this force, the impurity drifts and 
energy is dissipated. We assume that the impurities 
achieve their terminal drift velocity on a timescale 
much shorter than the frequencies of interest, and use 
Einstein's relation 112 



Dfi = Vik B T 
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to compute the drift velocity [T3]. Here D is the diffusion 
coefficient. The dissipated power per particle is < f^Vi >; 



summing over the particles and substituting the result 
into Eq. Q yields 
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Taking values of Suprasil 311 SV used in GEO-600 for 
OH ions [13], P(OH) = 3.9 • 10 24 m~ 3 and a = -4.52 • 
10~ 31 m 3 and estimating a value for D from [15] to be in 
the order of 10 _20 m 2 /s at room temperature we get 



S h (u) = 4.4 



10 _26 Hz 



-1/2 



(15) 



at a frequency of 100 Hz, suggesting that a simple version 
of the thermo-chemical noise is outside the realm of rele- 
vance for the GEO600 interferometer. However, there 
may be other optically-active mobile impurities (e.g., 
small structural defects or localized 2-state systems) in 
glass that have not yet been considered. Potential pres- 
ence of such impurities must be thoroughly investigated 
in future work. 



IV. DISCUSSION 

The calculations presented here demonstrate that the 
thermo-refractive fluctuations result in a noise floor 
which, while substantially below currently measured 
noise, may become important for future updates of the 
GEO600 interferometer. The standing-wave contribution 
should be taken into account in all future calculations of 
the thermo-refractive noise from transmissive optics. 

As a side product, we have identified a new type of 
thermal noise: the thermo-chemical noise, which is also 
enabled by the presence of the optical standing wave in 
the beamsplitter. The naive estimates of this noise, which 
are based on what is known about the chemical impu- 
rities in Suprasil 311, place it beyond the realm of con- 
cern. However, one needs to be vigilant about other types 
of mobile impurities which have not yet been identified, 
and which could potentially produce a higher thermo- 
chemical noise. 
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FIG. 2: Plot of measured GEO600 noise (black line), beamsplitter thermo-refractive noise from previous calculations [4] (red 
line) and corrected thermo-refractive noise using the same parameters (green line), and finally of corrected thermo-refractive 
noise using the updated parameters from [6] (blue line, [7]) 



